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p • Abstract 

' The matrix elements of electroweak currents which occur in exclusive decays of heavy 

hadrons are evaluated in the nonperturbative light-cone Fock representration. In 
general, each semileptonic exclusive decay amplitude receives two contributions, a 
diagonal An = parton-number-conserving amplitude and a An = —2 contribution 
in which a quark and an antiquark from the initial hadron Fock state annihilate to 
the leptonic current. The general formalism can be used as a basis for systematic 
approximations to heavy hadron decay amplitudes such as hard perturbative QCD 
contributions. We illustrate the general formalism using a simple perturbative model 
of composite hadrons. Our analysis demonstrates the occurence of "zero-mode" end- 
point contributions to matrix elements of the "bad" j~ currents in the Drell-Yan 
frame when —>■ 0. 

(Submitted to Physics Letters B.) 
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1 Introduction 



One of the most challenging problems at the intersection of quantum chromodynam- 
ics and electroweak physics is the evaluation of exclusive decay amplitudes of heavy 
hadrons such as the semileptonic decay n^i^u. The physics of such heavy 

hadron electroweak decays involve operator matrix elements which depend in detail 
on the quark and gluon composition of the initial state and final state hadrons. Even 
the presence of a heavy quark in the initial and/or final state does not simplify the 
complexity of the QCD analysis, since we must deal generally with hadron wavefunc- 
tions describing an arbitrary number of quark and gluon quanta. 

In this paper we shall give formulas for the current matrix elements {A\J^\B) 
describing general transition between hadrons B and A. The formulas are in prin- 
ciple exact, given the light-cone wavefunctions of hadrons. Our results generalize 
the expressions for the elastic form factors obtained by Drell and Yan 0, ^ and 
West 10]. The underlying formalism is the light-cone Hamiltonian Fock expansion in 
which hadron wavefunctions are decomposed on the free Fock basis of QCD. In this 
formalism, the full Heisenberg current can be equated to the current of the 
non-interacting theory which in turn has simple matrix elements on the free Fock 
basis. In the case of one-space and one-time theories, such as coUinear QCD [Q, the 
complete hadronic spectrum and the respective Fock state expansion can be deter- 
mined, at least numerically, using the DLCQ (Discretized Light-Cone Quantization) 
method [|5|. Eventually full solutions can be envisaged for physical theories such as 
QCD(3+1) using DLCQ, Wilson's front-form formalism, lattice analyses, and other 
non-perturbative Hamiltonian methods. For a review see Ref. 0. 

An exact formalism provides the opportunity to make systematic approximations 
and account for negelected terms. For example, we can identify the contributions to 
exclusive decay amplitude of heavy hadrons which can be accounted for by hard per- 
turbative QCD effects On the other hand, we also can identify specific physical 
mechanisms which are due to the presence of higher Fock state non-valence configu- 
rations of the hadrons. 

It is well known that elastic form factors at space-like momentum transfer 
= —Q^ < are most simply evaluated from matrix elements of the "good" current 
j+ = + in the preferred Lorentz frame where g"*" = + = 0. The current 
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has the advantage that it does not have large matrix elements to pair fluctuations, so 
that only diagonal, parton-number-conserving transitions need to be considered. The 
use of the j"*" current and the q'^ — frame brings out striking advantage of the light- 
cone quantization formalism: only diagonal, parton-number-conserving Fock state 
matrix elements are required. However, in the case of the time-like form factors 
which occur in semileptonic heavy hadron decays, we need to choose a frame with 
> 0, where q'^ is the four-momentum of the lepton pair. Furthermore, in order 
to sort out the contributions to the various weak decay form factors, we need to 
evaluate the "bad" — current j~ = j° — as well as the "good" current j^. In such 
cases we will also require off-diagonal Fock state transitions; i.e. the convolution of 
Fock state wavefunctions differing by two quanta, a qq' pair. The entire electroweak 
current matrix element is then in general given by the sum of the diagonal n ^ n 
and off-diagonal n + 1 ^ n — 1 transitions. As we shall see, an important feature of 
a general analysis is the emergence of singular 6{x) "zero-mode" contributions from 
the off-diagonal matrix elements if the choice of frame dictates g"*" — 0. 



2 Matrix Elements of Electroweak Currents 

The light-cone Fock expansion is defined as the projection of an exact eigensolution 
of the full light-cone quantized Hamiltonian on the solutions of the free Hamiltonian 
with the same global quantum numbers. The coefficients of the Fock expansion are the 
complete set of n-particle light-cone wavefunctions, {ipnixi, k±i, Xi)}. The coordinates 
Xi, k±i are internal relative coordinates, independent of the total momentum of the 
bound state, and satisfy < Xi < 1, J2i = 1 and = 0^- Here x — — 

po'^p3 and we use the metric convention a ■ b = ^{a'^b~ + a~b'^) — a± ■b±. 

The evaluation of the semileptonic decay amplitude B AiV requires the matrix 
element of the weak current between hadron states (i?|j^(0)|A). (See Fig. 1.) The 
interaction current then has simple matrix elements of the free Fock amplitudes, with 
the provisal that all Xi > 0. We shall adopt the choice of a Lorentz general frame 
where the outgoing leptonic current carries = {q^, q±, q~) = ^AP+, q±, '^^+ ^ . In 
the limit A — 0, the matrix element for the -|- vector current should coincide with 
the Drell-Yan formula. 
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Figure 1: Exact representation of electroweak decays and time- like form factors in 
the light-cone Fock representation. 

For the n — > n diagonal term (An = 0), the final-state hadron wavefunction 
has arguments , k±i — jz^q± for the struck quark and ^z^, k±i + ■^z^q± for the 
n—1 spectators. We thus have a formula for the diagonal (parton-number-conserving) 
matrix element of the form: 

Mk-iBu.o = i:n/>-.i'^-.<..i./^*(i-|:-.) ^-'''(eSx,^ 



where 



x' 



A 1 

k'j_i — k±i — jza^-^ struck quark 



l-A ' 



(2) 



, k'j_^ = k±i + ■jz^q± for the (n — 1) spectators. 
A sum over all possible helicities Aj is understood. If quark masses are neglected 
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the vector and axial currents conserve helicity. We also can check that I^r^i — 

For the n + 1 — > n — 1 off-diagonal term (An = —2), let us consider the case where 
partons 1 and n + 1 of the initial wavefunction annihilate into the leptonic current 
leaving n — 1 spectators. Then x„+i = A — Xi, k±n^i = q± — k±i. The remaining 
n — 1 partons have total momentum ((1 — A)P+, —q±). The final wavefunction then 
has arguments x'^ = (YZa) ^J-* ~ ^-L* + T^^^- O^ViS obtain the formula for 
the off-diagonal matrix element: 
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T^fJ 2(27r)3 Ijio V 2(2 

n+l \ /n+1 

i=i / Vi=i 

><V'A(n-l)(^i'^li'^i) f i^B{n+l){{Xl,X,,Xn+l = A - xj, 

{fc^i, /c_L„+i = q± — k±i}, {Ai, Aj, A„_|_i = — Ai}). (3) 



Here i = 2, 3, ■ ■ ■ , n with 

label the n — 1 spectator partons which appear in the final-state hadron wavefunc- 
tion. We can again check that the arguments of the final-state wavefunction satisfy 

Y17=2 = I5 J27=2 ^'±i — 0^- 

The free current matrix elements j'^ in the light-cone representation are easily 
constructed. For example, the vector current of quarks takes the form 

^ u{x',k'^,X')-ff'u{x,k^,X) 

^/k+^/k+' 

and 

J+ = 25a,a' . 

The other light-cone spinor matrix elements of j'^ can be obtained from the tables in 
ref. [|]. In the case of spin zero partons 

. , X + x' 
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and 

_k- + k'- 

However, instead of evaluating each k~ in the j~ current from the on-shell condition 
k~k~^ = m?, one must instead evaluate the k~ of the struck partons from energy 
conservation k~ = Pj^itiai ~Pspectator- This effect is seen explicitly when one integrates 
the covariant current over the denominator poles in the k~ variable. It can also be 
understood as due to the implicit inclusion of local instantaneous exchange contribu- 
tions obtained in light-cone quantization j^, ITO). The mass rn'^pccta.toT which is needed 



for the evaluation of j current in the diagonal case is the mass of the entire specta- 
tor system. Thus ml^^^^^^^^^ = m^^^^^^^^^ + A^ispectator, where fc±spectator = Ej k±j and 
"^ispectator /^spectator = J^j^'j/^j, summcd ovcr the j spcctators. This is an impor- 
tant simplification for phenomenology, since we can change variables to ?T^spectator ^^^1 
o?^fc±spectator and replace all of the spectators by a spectral integral over the cluster 
mass /^spectator- ^ specific example is presented in the next section. 



3 Example — 0^ Perturbation Theory 

As an explicit example and check on the above formalism, we shall consider the elec- 
tromagnetic vector current matrix element of a neutral composite system composed 
of two charged scalars where the light-cone wavefunctions are known explicitly from 
perturbation theory. To construct the model, we consider a 3+1 dimensional system 
represented by the Lagrangian: 

C = {d^(f)a + ieaAMHd''<l)a + ieaA''(f)a) -ml(j)l(j)a (5) 
+ {df,(j)b - iebA^(j)i,)\d^(f)b - iebA^b) - ml4>l(Pb 

The composite system wavefunction can be normalized to unity by a choice of the 
effective coupling g. 

We can derive the light-cone amplitudes from the covariant amplitude by inte- 
grating over the k~ variable The amplitude of the process drawn in Fig. 2 is 
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Figure 2: Scalar perturbation theory model for semileptonic decay, (a) Covariant 
representation, (b), (c) Light-cone time-ordered contributions to the decay ampli- 
tude. These contributions can be identified as the convolution of light-cone Fock 
wavef unctions with 2^2 and 3^1 parton number, respectively. 



given as follows from the Feynman rules: 



Mt 



leag 



(6) 



{2k - qY 



(27r)4 (A;2 - + ie) {{k - qf -ml + ie){{k - Py - ml + ie) 
cPk^ /• „^ , 1 



2(27r)^ 



P+dx 



X / dk~ 



{2k - qY 



{k- - q-) - 



(mg + (fcx-g'±)^-ie 
ix-A)P+ 



{k--P-) + 



(mg+fcx)^)-^6 ^ 
{l-x)P+ ^ 



where we used k^ = xP^ . When we perform the integration over k , the integral 
does not vanish only for < x < 1. 
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For A < X < 1, which corresponds to Fig. 2(b), 

{Mt,M-,M^a)i2-.2) (7) 



2(2 



TT 



X I (2x - A) , ^ I 2 I - 1 - 1 . (2., - 



1 

X 



x{x — A)(l — a; 

X 



a; J y 1-A x-A 1-x 



1 r ^2 



2(27r)2 



x|(2.-A),^(2(M^-M±|i)l-(il±ii)|,(2..-f.) 

(1-A) 

X- 



■a;(x - A)(l - x) 

X 



P+e. 



ag^ J^dx J 



X 1—x 



1-A 1-A 



2(27r) 



3 



1111 

X 



A 



_ — * 



2(27r)3 

X'0(2)«, 4' ^'a±: k'b±; M, m„, m6)jf2^2)aV'{2)(2;a, Xb, Kx, hi.] M, Ula, Ulb) , 



where 



J(2^2)a = ea , (8) 

— Aj 

X I (2. - A, . ^ f2 f - ^±^] - nil I , p.-, 
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-0(2) {Xa, Xb, ka±, kb±] M, nia.nib) = Q 



\ Xa Xt, J 

Xa = X, Xfc = 1 X, ka± = k±, kb± = k±, 

X'a = ^ , ^b^Y~r^' k'a±^ka± + (l-x'JP'±, k'b± ^ kb± - x'^P' ±, 

P'± = -q±. 

For < X < A, which corresponds to Fig. 2(c), 

{M^,M-,M^a)i3^1) (9) 



dpk\ 



2(27r)3 

x|(2.-A),^f2M±^i)-(^| ,(2*.-fJ 

1 



X 



X 



x(A -x){l- x) 

X 



I X 1—x J \ X A—x 1— A J 



xtp{i) L _ ^ , ky± - P'± ] j^s^i)a'^{3){x, y,A-x, k±, ky±, q± - k±; M, rUa, rUb), 
where P'± — —q± and 



" ^x{^-x) 



x((2.-A),^(2(!^-(^| ,(2.^,-,l)|, (10) 

^/'(3)(a;, y,A-x, k±, ky±, q_L - k±; M, nia, rrib) 
= 9 



1 



x{l — x)^(A — x)y 



X 



X 1—x I \. X A—x y I 



M2 



V'(i)(i/,M = ^^<^(y-i)5'(M- 
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Then, by combining (|^) and (H) the amphtude is given as: 



1 r 

X {Xa, Xh, ka±, h±; M, rUa, 171^) 

+ dx TTTT— / dy 



J 2(271)3 70 2(27r)3 
x^(i) V - j(3^i)a'>Pi3)ix, y,^-x,k±, ky^, q±-kr, M, rUa, mb). 

By adding the above amphtude Ai^ and that given by exchanging a and b {ca + 
Cft = 0), we obtain the total amphtude: 

= {2P-qfF{q^) (12) 

= J^dx J ^^^^^3 ^(2)«, a;;, A;'a±, M, rUa, "^6)j('2^2)a^(2) 

X (xa, Xb, ka±, kb±; M, rria, mb) 



J 2(27r)3 Jo 2(27r)3 

X^(l) (y^A ' ~ ^'^) ^r3-l)a^(3) (a;, A - X, fcx, <?1 - ^, "^a, "^b) 



+ (a ^ 6 



where = (2P - qYF{q^) foUows from q^M^ = 0. 

For q^ ^ 0, q± ^ 0± and A ^ 0, + component of (|12D gives 



/(Z njj_ f -* -* 

— --^ |V'(2)(a;, 1 - x, fc_L, M,ma,mb)p (13) 



(27r)3 Jo 

(ix |^/'(2)(x, 1 -X, fc_L,-fc_L;M,m6,ma)|^ = 0, 



where 'ip(2) is given in (H). Each term can be normahzed to unit charge, thus providing 
wavefunction renormahzation in the modeL Ahernatively we can evaluate the — 



component of ([T^) to obtain 
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m (14) 



d'k^ 1 / /•! , 1 ml + kl 



1 2^ 



H 



5 



where the J^"^^ and -^t^q"^ terms come from the singular contributions of the 
/q^ ^^(3) terms in (|12[) when we take the limit A ^ 0. The ± components 

of (|T^) do not give more information since (2P — ^± — > 0_l in the left hand side and 
the integrand of the right hand side is odd about k±. 

The above analysis provides an explicit realization of the general formulas (|l|) and 
(13). In this simple model two transition matrix elements appear: 2 — 2 and 3-^1. 



The equality of the formulas for (^) and ([1^ is a general condition which follows 
from gauge invariance and consistency of the light-cone formalism. We have verified 
the equality for the perturbative model by direct evaluation of the integrals. 

In the case of general composite systems, the equality of the form factors at zero 
momentum transfer obtained from the and J~ currents provides a type of virial 
theorem for the matrix elements ([1|) and (|^). In general the two determinations of 
the total charge F{q^ = 0) must be consistent: 

F(0) = -4^ lim (a\J+\b) , (15) 
F(0) = lim ((A\J-\b) +(A\J-\b) ). (16) 

^ ^ 2P~ q^q^A^O \\ ' ' /An=0 \ ' ' / An=-2 J ^ ^ 

Here P^P~ = M^. Note that the second term of (^) includes the zero mode 6{x) 
contributions from the n + 1 — n — 1 off-diagonal matrix element. 
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4 Conclusions 



A most important feature of the light-cone formahsm is that all matrix elements 
of local operators can be written explicitly in terms of simple convolutions of light- 
cone Fock wavefunctions {ipni^i, kj_i, Xi)}. In the case of exclusive semileptonic B- 
decays, such as B ^ Aiu, the decay matrix elements require the computation of 
the diagonal matrix element n n where parton number is conserved and the off- 
diagonal n + 1 — s> n — 1 convolution where the current operator annihilates a qq' 
pair in the initial B wavefunction. This term is a consequence of the fact that the 
time-like decay = {pe + Prr)'^ > requires a positive light-cone momentum fraction 
> 0. Conversely for space-like currents, one can choose g"*" = 0, as in the Drell- 
Yan-West representation of the space-like electromagnetic form factors. However, as 
we have seen from the explicit analysis of the form factor in a perturbation model, 
the off-diagonal convolution can yield a nonzero q^/q~^ limiting form as q'^ —>■ 0. This 
extra term appears specifically in the case of "bad" currents such as J~ in which 
the coupling to qq fluctuations in the light-cone wavefunctions are favored. In effect, 
the q^ —* limit generates 6{x) contributions as residues of the n + 1 n ~ 1 
contributions. The necessity for this zero mode 6{x) terms were first noted in the 
pioneering work of Chang, Root and Yan |p, and Burkardt analyzed it in his studies 
of higher-twist parton distributions []T2|. Here we see that the presence of such terms 



are a general feature of local operator matrix elements when one selects the simplified 
g"*" = frame. 

We have also seen that the proper treatment of the J~ current implies new consis- 
tency conditions which must be obeyed by the light-cone wavefunctions. For example, 
current conservation for the form factors of spin zero hadrons requires 

i2p-qrF{q') = {p-q\riO)\p) (17) 

and thus 

{P-^\J^\P)=^^^ {P-^\J'\P) ■ (18) 

We have explicitly verified this new type of virial theorem in a simple scalar composite 
model in section 3. 

The off-diagonal + 1 — > — 1 contributions provide a new perspective on the 
physics of i?-decays. A semileptonic decay involves not only matrix element where 
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a quark changes flavor, but also a contribution where the leptonic pair is created 
from the annihilation of a qq' pair within the Fock states of the initial B wavefunc- 
tion. The semileptonic decay thus can occur from the annihilation of a nonvalence 
quark-antiquark pair in the initial hadron. This feature will carry over to exclusive 
hadronic -B-decays, such as n^D^. In this case the pion can be produced from 

the coalescence of a du pair emerging from the initial higher particle number Fock 
wavefunction of the B. The D meson is then formed from the remaining quarks after 
the internal exchange of a 1^ boson. 

We have emphasized the remarkable advantage of the light-cone formalism that 
all matrix elements of local operators can be written down exactly in terms of simple 
convolutions of light-cone Fock wavefunctions. The light-cone wavefunctions depend 
only on the hadron itself; they are pro cess- independent. The formalism is relativis- 
tic and frame-independent — the incident four-vectors can be chosen in any frame. 
Note that the matrix element of a current in the covariant Bethe-Salpeter formalism 
requires the construction of the current from insertions into an infinite number of 
irreducible kernels. The Bethe-Salpeter formalism becomes even more intractable for 
bound-states of more than two particles. 

In principle, a precise evaluation of the hadronic matrix elements needed for B- 
decays and other exclusive electroweak decay amplitudes requires knowledge of all of 
the light-cone Fock wavefunctions of the initial and final state hadrons. In the case 
of some model gauge theories such as QCD or coUinear QCD in one-space 
and one-time dimensions, the complete evaluation of the light-cone wavefunction is 
possible for each baryon or meson bound-state using the discretized light-cone quan- 
tization method. It would be interesting to use such solutions as a model for physical 
5-decays. 

The evaluation of the light-cone Fock wavefunctions in QCD (3+1) is not at present 
computationally feasible because of the large number of degrees of freedom within the 
hadron wavefunctions. Nevertheless, the existence of an exact formalism provides a 
basis for systematic approximations and a control over neglected terms. For example, 
one can analyze exclusive semileptonic 5-decays which involve a heavy internal mo- 
mentum transfer using a perturbative QCD formalism patterned after the analysis of 
form factors at large momentum transfer |^. The hard-scattering analysis proceeds 
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by writing each hadronic wavefunction as a sum of soft and hard contributions 



^Pn = rf\Ml < A') + ^P'^^'^Ml > A'), (19) 

where 

^HF^t ,20) 

is the invariant mass of the partons in the n-particle Fock state and A is the separation 
scale. The high internal momentum contributions to the wavefunction ipn'^'^'^ can be 
calculated systematically from QCD perturbation theory from the interaction of the 
gluon exchange kernel. The contributions from high momentum transfer exchange 
to the -B-decay amplitude can then be written as a convolution of a hard scattering 
quark-gluon scattering amplitude Th with the distribution amplitudes (j){xi,A), the 
valence wavefunctions obtained by integrating the constituent momenta up to the 
separation scale A4n < A < Q. This is the basis for the perturbative hard scattering 
analyses of Refs. |14|, |1^, [1^. In our exact analysis, one can identify the hard PQCD 



contribution as well as the soft contribution from the convolution of the light-cone 
wavefunctions. Furthermore, the hard scattering contribution can be systematically 
improved. For example, off-shell effects can be retained in the evaluation of Th by 
utilizing the exact light-cone energy denominators. This effect will be analyzed in a 
separate paper. 
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